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Abstract 

Observations are reported in isothermal torsional oscillation tests on melts of isotactic 
polypropylene (iPP) and low-density polyethylene (LDPE) in the intervals of tempera- 
ture between 190 and 250 °C (iPP) and between 120 and 190 °C (LDPE). With reference 
to the concept of transient networks, constitutive equations are developed for the vis- 
coelastic response of polymer melts at three-dimensional deformations with small strains. 
A melt is treated as an equivalent network of strands bridged by temporary junctions 
(entanglements and physical cross-links whose life-times exceed the characteristic time of 
deformation). The time-dependent behavior of the network is modelled as detachment 
of active strands from their junctions and merging of dangling strands with the network. 
The network is assumed to be strongly heterogeneous in the sense that different junctions 
have different activation energies for separation of strands. The stress-strain relations 
involve three adjustable parameters (the plateau modulus, the average activation energy 
for rearrangement of strands and the standard deviation of activation energies) that are 
determined by matching the dependencies of storage and loss moduli on frequency of os- 
cillations. The difference in the effects of temperature on the material constants of iPP 
and LDPE is associated with the difference in their molecular architecture. 
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1 Introduction 



This paper is concerned with the effect of temperature on the viscoelastic response of melts 
of isotactic polypropylene (iPP) and low-density polyethylene (LDPE). The choice of these 
polyolefins for the investigation may be explained by two reasons. First, polyethylene and 
polypropylene are conventional polymers widely used in industrial applications: oriented films 
for packaging, reinforcing fibres, non- woven fabrics, pipes, etc., as well as components of copoly- 
mers and blends with improved mechanical properties. Secondly, these polymers have a similar 
structure of backbones (polypropylene differs from polyethylene by the presence of methyl side- 
groups only), but quite different molecular architecture (relatively short-branched chains in iPP 
versus long-branched chains in LDPE). The aim of the present work is to shed some light on re- 
lations between the micro-structure of an ensemble of macromolecules on the one hand, and the 
time-dependent behavior of polymer melts observed in conventional torsional oscillatory tests, 
on the other. In particular, we focus on the effects of long chain branches on the evolution of 
(i) the shear modulus of a melt and (ii) its characteristic relaxation time with temperature. 

According to the theory of entropic elasticity (Ferry, 1980), the elastic modulus p of a 
flexible chain linearly increases with the absolute temperature T, 

p = p Q k B T, (1) 

where k B is Boltzmann's constant, and the dimensionless coefficient p is of order of unity (this 
coefficient is determined by an averaging method in the calculation of the chain's entropy). As 
conventional concepts in rubber elasticity disregard interactions between chains (these inter- 
actions are accounted for by means of the incompressibility condition), the shear modulus of 
an ensemble of polymer chains equals the product of the concentration of strands (segments of 
macromolecules between contiguous junctions) per unit mass N by the elastic modulus of an 
individual chain p, 

G = ppN, (2) 

where p stands for mass density. For linear chains, the number of strands per unit mass iV 
equals the product of the number of chains per unit mass Af by the average number of strands 
per chain f3, 

N = (3Af , j3=±-. (3) 

Here L is the average contour length of a chain, and L e is the average contour length between 
entanglements. Combining Eqs. (1) to (3), we arrive at the conventional formula 

G = pp k B AfL^. (4) 

Experimental data demonstrate a pronounced decrease of the apparent elastic modulus (de- 
fined, e.g., as the storage modulus G' measured at the maximal frequency cj max available in an 
experiment) with temperature T (see Figures 1 and 3 below). As all parameters on the right- 
hand side of Eq. (4) but the ratio T/L c are temperature independent (insignificant changes in 
density with temperature are disregarded in the present study), this implies that the average 
contour length between entanglements L e strongly increases with temperature, 

L e (T) = Cl c (T), Z e (T) = -^, C = pp Q k B ML. (5) 
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For an ensemble of linear chains, the growth of L e with temperature T may be attributed to an 
increase in the average inter-chain distance driven by activation of thermal fluctuations. The 
same process determines an increase in the specific volume (length) with temperature, which 
is described by the conventional equation 

IdT = (6) 
where / denotes a characteristic length, and a is the coefficient of linear thermal expansion. 
Based on this similarity, we propose to describe changes in l e {T) by the differential equation 
similar to Eq. (6), 

i dl e ^ 

l e -lf)dT~ a ^ [n 

where a e stands for an analog of the coefficient of linear thermal expansion for the average 
length between entanglements, and characterizes the contour length between entanglements 
at relatively low temperatures T. Solving Eq. (7) for a constant coefficient a c , we find that 

/e(T)=/f + /«exp(^), T e = l. (8) 



Combining Eqs. (5) and (8), we arrive at the formula for the shear modulus G as a function of 
temperature T, 

G(T) = /(°) + /W T (t/tV (9) 
le ' + le exp(T/T c ) 

It is worth noting that an equation similar to Eq. (9) (where the coefficient l^ is absent) has 
recently been proposed by Wood-Adams and Costeux (2001) based on another physical ground. 

A substantial increase in the entanglement distance l e with temperature T has been observed 
by Richter et al. (1990, 1993) about a decade ago by using neutron spin-echo spectroscopy 
(NSES) in polyisoprene, polybutadiene and a poly(ethylene-propylene) alternative copolymer. 
To the best of our knowledge, no attempts have been made, however, to model the growth of 
Z e (T) with the help of differential equations, as well as to evaluate their parameters by matching 
observations in rheological tests. The latter may be explained by the use of conventional 
methods (Ferry, 1980) for determining the plateau modulus G° (an analog of the instantaneous 
modulus G). These methods are grounded on shifting the graphs of the storage modulus 
G' measured at various temperatures T and plotted in the double-logarithmic coordinates as 
functions of frequency u along both axes: no restrictions are imposed on the shift along the 
abscissa axis (which characterizes the shift factor), whereas an appropriate shift along the 
ordinate axis is presumed to be proportional to the logarithm of temperature T, which is 
tantamount to formula (2) with a constant N. 

The objective of this study is two-fold: 

• to develop a constitutive model for the isothermal time-dependent response of a polymer 
melt at small strains that involves only three material constants (including the instanta- 
neous shear modulus G), 

• to report experimental data in torsional oscillation tests on isotactic polypropylene and 
low-density polyethylene at various temperatures T, to determine adjustable parame- 
ters in the stress-strain relations by fitting the observations, and to verify Eq. (9) by 
comparing the experimental data with the results of numerical simulation. 
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To develop constitutive equations tractable from the mathematical standpoint, we adopt a 
homogenization concept. According to it, a complicated micro-structure of a polymer melt is 
replaced by an equivalent phase, whose response captures essential features of the mechanical 
behavior. Following common practice (Sweeney et al., 1999), a network of chains is chosen as 
an equivalent phase. 

The time-dependent behavior of melts is modelled within the concept of transient networks 
(Green and Tobolsky, 1946; Yamamoto, 1956; Lodge, 1968; Tanaka and Edwards, 1992). A 
polymer melt is thought of as a network of strands bridged by temporary junctions (entan- 
glements and physical cross-links whose life-time exceeds the characteristic time of rheological 
tests). It is assumed that active strands (whose ends are linked to contiguous junctions) sepa- 
rate from these junctions, while dangling strands (that have a free end) merge with the network. 
Detachment and attachment events occur at random times as appropriate strands are excited 
by thermal fluctuations. Following Drozdov and Christiansen (2003), we assume the network 
to be strongly heterogeneous in the sense that different junctions have different activation en- 
ergies for detachment of active strands. The inhomogeneity of the network is attributed to 
(i) density fluctuations in the ensemble of macromolecules, and (ii) the fact that the life-time 
of entanglements created by long branches as well as formed by backbones but located in the 
neighborhoods of chain's ends is smaller than that of knots between backbones of different 
macromolecules entangled in the vicinity of their middle-points. The distribution of active 
strands with various activation energies for rearrangement is described by the random energy 
model (Derrida, 1980). 

The stress-strain relations involve three adjustable parameters: (i) the instantaneous mod- 
ulus G, (ii) the average activation energy for separation of strands from temporary junctions 
V*, and (iii) the standard deviation of activation energies E*, which, in general, depend on 
temperature T. These quantities are found by fitting experimental data for the storage and 
loss moduli in shear oscillatory tests at each temperature separately. This procedure provides 
an opportunity to evaluate the dependence G(T) and to verify Eq. (9). 

The exposition is organized as follows. Experimental data in torsional oscillation tests 
are reported in Section 2. Constitutive equations for an heterogeneous transient network of 
strands at three-dimensional deformations are derived in Section 3. Adjustable parameters in 
the stress-strain relations are determined in Section 4. A brief discussion of our findings is 
given in Section 5. Some concluding remarks are formulated in Section 6. 

2 Experimental procedure 

Isotactic polypropylene PP 1012 (density 0.906 g/cm 3 , melt flow rate 1.2 g/10 min) was pur- 
chased from BP Amoco Polymers, Inc. Low-density polyethylene Huntsman PE 1020 (density 
0.923 g/cm 3 , melt flow rate 2.0 g/10 min) was supplied by GE Company. Granules were dried 
at the temperature T = 100 °C for 12 h prior to molding. Circular plates with radius 62 
mm and thickness 3 mm were molded in injection-molding machine Battenfeld 1000/315 CDC 
(Battenfeld). Specimens for rheological tests (with diameter 30 mm) were cut from the plates. 

To evaluate melting temperatures of iPP and LDPE, DSC (differential scanning calorime- 
try) measurements were carried out by using DSC 910S apparatus (TA Instruments). The 
calorimeter was calibrated with indium as a standard. Two specimens of each polymer with 
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weights of approximately 15 mg were tested with a heating rate of 10 K/min from room tem- 
perature to 200 °C. The melting temperatures T m = 172 (iPP) and T m = 113 °C (LDPE) were 
determined as the point corresponding to the peaks on the melting curves. 

Rheological tests were performed by using RMS-800 rheometric mechanical spectrometer 
with parallel disks (diameter 25 mm, gap length 2 mm) at the temperatures T = 190, 210, 
230 and 250 °C (iPP) and T = 120, 130, 140, 150, 160, 170, 180 and 190 °C (LDPE). Given 
a temperature T, at least two dynamic tests were carried out on different samples. The shear 
storage modulus G' and the shear loss modulus G" were measured in oscillation tests (the 
frequency-sweep mode) with the amplitude of 15 % and various frequencies u ranging from 0.1 to 
100 rad/s. Our choice of the amplitude of oscillations was based by the following requirements: 
(i) mechanical tests were performed in the region of linear viscoelasticity, and (ii) the torque was 
less than its ultimate value 0.2 N-m. The limitation on the minimum frequency of oscillations 
was imposed by two conditions: (i) the torque exceeded its minimum value 2.0 • 10~ 4 N-m, 
and (ii) the duration of a test did not exceed 20 min, which ensured that thermal degradation 
of polymers at elevated temperatures may be disregarded. To check that the storage and loss 
moduli were not affected by the strain amplitude, several tests were repeated with the amplitude 
of 5 %; no changes in dynamic moduli were observed. The temperature in the chamber was 
controlled with a standard thermocouple that indicated that the temperature of specimens 
remained practically constant (with the accuracy of ±1.0 °C). 

Each test was performed on a new sample. The specimen was thermally equilibrated in 
the spectrometer (during 5 min), the gap length was reduced to 2 mm, an extraneous material 
was carefully removed, and the storage and loss moduli were measured at various frequencies 
uj starting from the lowest one. 

The storage G' and loss G" moduli are depicted versus the logarithm (log = log 10 ) of 
frequency uj in Figures 1 and 2 for iPP and in Figures 3 and 4 for LDPE. Conventional semi- 
logarithmic plots are used to characterize changes in these quantities with frequency. According 
to Figures 1 to 4, the dependencies of storage and loss moduli on frequency of oscillations have 
similar shapes for both polyolefins. Given a temperature T, the storage modulus G' and the 
loss modulus G" strongly increase with frequency uj. For a fixed frequency uj, the dynamic 
moduli pronouncedly decrease with temperature T. 

3 Constitutive equations 

With reference to the concept of transient networks, a polymer melt is thought of as an equiv- 
alent network of strands bridged by temporary junctions. A strand whose ends are linked to 
contiguous junctions is treated as an active one. When an end of an active strand separates 
from a junction, the strand is transformed into the dangling state. When a free end of a dan- 
gling strand captures a nearby junction, the strand returns into the active state. Separation of 
active strands from their junctions and merging of dangling strands with the network occur at 
random times when the strands are excited by thermal fluctuations. According to the theory of 
thermally-activated processes (Eyring, 1936), the rate of detachment of strands from temporary 
junctions T is governed by the equation 




(10) 
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where T is the attempt rate (the number of separation events per strand per unit time), k-Q 
is Boltzmann's constant, T is the absolute temperature, and v > is the activation energy for 
separation of an active strand. In what follows, we set T = 10 11 s _1 , which corresponds to the 
characteristic relaxation rate at the monomeric scale (deGennes, 1979). 

For isothermal deformation at a temperature T, we introduce the dimensionless activation 
energy 

and present Eq. (10) in the form 

r» =r exp(-t;). (12) 

To describe the time-dependent response of a melt, we suppose that different junctions are 
characterized by different dimensionless activation energies v (Drozdov and Christiansen, 2003). 
The distribution of active strands in a transient network is determined by the number of active 
strands per unit mass N and the distribution function p(v). The quantity Np{v)dv equals the 
number of active strands per unit mass linked to junctions with the dimensionless activation 
energies u belonging to the interval [v, v + dv]. 

Separation of active strands from temporary junctions and merging of dangling strands with 
the network are entirely described by the function n(t, r, v) that equals the number (per unit 
mass) of active strands at time t > linked to temporary junctions with activation energy v 
which have last merged with the network before instant r e [0,t]. 

The quantity n(t, t, v) equals the number of active strands (per unit mass) with the activa- 
tion energy v at time t, 

n(t,t,v) = Np(v). (13) 



The function 



(14) 



determines the rate of reformation for dangling chains: the amount 7(r, v)dr equals the number 
of dangling strands (per unit mass) that merge with temporary junctions with activation energy 
v within the interval [r, r + dr] . The quantity 

dn . . 
— (t,T,v) dr 

is the number of these strands that have not separated from their junctions during the interval 
[r, t] . The amount 

-—(t,0,v) dt 

is the number of active strands (per unit mass) that detach (for the first time) from the network 
within the interval [t,t + dt], while the quantity 

■(t, r, v) dtdr 



dtdr 



equals the number of strands (per unit mass) that have last merged with the network within 
the interval [r, r + dr] and separate from the network (for the first time after merging) during 
the interval [t, t + dt]. 
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The rate of detachment T is defined as the ratio of the number of active strands that separate 
from temporary junctions per unit time to the total number of active strands. Applying this 
definition to active strands that were connected with the network at the initial instant t — 0, 
and to those that merged with the network within the interval [r, r + dr\ , we arrive at the 
differential equations 



dn 
~dt 



(t,0,v) = -F(v)n(t,0,v), 



d 2 n , . . s dn . . 

- (t ,r,„) = -r W -( t ,r,„). 



(15) 



Integration of Eq. (15) with initial conditions (13) (where we set t — 0) and (14) implies that 

dn 



n(t,0,v) = Np(v)exp —T(v)t 



dr 



(t, t, v) = 7(r, v) exp —T(v) (t — t) 



To exclude the function 7(t, v) from Eq. (16), we use the identity 

rt Qji 

n(t,t,v) — n(t,0,v) + —(t,r,v)dr. 

Jo or 

Substitution of expressions (13) and (16) into Eq. (17) results in 

Np(v) = Np(v)exp\-r(v)i\+ I t(t, v) exp\-F(v)(t - t) 
L J jo L 



dr. 



(16) 



(17) 



(18) 



The solution of linear integral equation (18) reads j(t,v) = Np(v)T(v). It follows from this 
equality and Eq. (16) that 



dn 

— {t, r, v) = Np(v)r(v) exp[-I»(* - r) 



(19) 



We adopt the conventional assumptions that (i) the excluded-volume effect and other multi- 
chain effects are screened for individual strands by surrounding macromolecules, (ii) the energy 
of interaction between strands can be taken into account with the help of the incompress- 
ibility condition, and (iii) thermal oscillations of junctions can be disregarded, and the strain 
tensor for the motion of junctions at the micro- level coincides with the strain tensor for macro- 
deformation. 

At isothermal deformation with small strains, a strand is treated as an isotropic incom- 
pressible medium with the strain energy 



w = lie : e , 

where e is the strain tensor for transition from the reference (stress-free) state of the strand to 
its deformed state, the average elastic modulus p, is given by Eq. (1), the prime stands for the 
deviatoric component of a tensor, and the colon denotes convolution of two tensors. 

According to the affinity hypothesis, the strain energy wo(t, 0) of an active strand that has 
not separated from the network during the interval [0, t] reads 

w(t,0) = pe{t) : e'(t), 

where e(t) is the strain tensor for transition from the initial (stress- free) state of the network 
to its deformed state at time t. With reference to Tanaka and Edwards (1992), we suppose 



7 



that stress in a dangling strand totally relaxes before this strand captures a new junction. This 
implies that the stress-free state of an active strand that merges with the network at time r > 
coincides with the deformed state of the network at that instant. The strain energy of an active 
strand that has last merged with the network at time r e [0, t] is given by 

w(t,r)=p[e(t)-e(r)]': [e(t) - e(r) ' 

Multiplying the strain energy per strand by the number of active strands per unit mass and 
summing the mechanical energies of active strands linked to temporary junctions with various 
activation energies, we find the strain energy per unit mass of an equivalent network 

fOO ( ft f)nr, r 1 / r 1 / "1 

W(t) = /I J |n(t, 0, v)e'(t) : e'(t) + J — (t, r, v) [e(t) - e(r)] : [e(f) - e(r)J drjdv. (20) 

Differentiating Eq. (20) with respect to time t and using Eqs. (16) and (19), we arrive at the 
formula 

™(t) = A'{t):^{t)-B{t), (21) 



where 



f ft fOO 

A(t) = 2fjiN^i(t)-j e(r)dr J T{v)exp 
B{t) = nj™r(v){n(t,0,v)?(t):?(t) 



-T(v)(t — t) p(v)dv >, 



ft i ' r i ' 1 

+ / -z-(t,T,v) e(t) -eir) : eit) - eir) dr\dv>0. 
Jo or 1 J L J J 



(22) 



(23) 



For isothermal deformation of an incompressible medium, the Clausius-Duhem inequality reads 

Q - _^ + ^ ■ ^ > 

dt p dt ~ ' 

where Q is internal dissipation per unit mass, and a stands for the stress tensor. Substitution 
of Eq. (21) into this equation implies that 



Q(t) = -[a'(t)-pA\t)]:^(t) + B(t)>0. 
p l 1 dt 



(24) 



As the function B(t) is non- negative, see Eq. (23), dissipation inequality (24) is satisfied, 
provided that the expression in the square brackets vanishes. This assertion together with Eq. 
(22) results in the constitutive equation 



a(t) = -P(t)I + 2G{e'(t) - e'(T)dT F(v) exp[-I»(* - r)]p(u)du|, 



(25) 



where P(t) is pressure, / is the unit tensor, and the shear modulus G is determined by Eq. (2). 

Formula (25) describes the time- dependent response of an equivalent network at arbitrary 
three-dimensional deformations with small strains. This equation implies that in a shear test 
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with e(t) = e(t)eie 2 , where e(t) is the shear strain, and e m (m = 1,2,3) are unit vectors of a 
Cartesian frame, the shear stress a(t) reads 



r rt roc f -| 

a(t) = 2GU(t) - J e{r)dT J F(v) exp[-I»(i - T)\p{v)dv 



(26) 



It follows from Eq. (26) that in a shear oscillation test with e(t) = e exp(iut), where e and uj 
are the amplitude and frequency of oscillations, and i = the transient complex modulus 

G*(t,u) = a(t)/(2e(t)) is determined by the formula 

G*(t,u) = cjl - J™ T{v)p{v)dv J\xp[-(r(v) +io;)s]ds|, 

where s = t — r. This equality implies that the steady-state complex modulus G*(cu) = 
lim^oo G*(t, uj) is given by 

G*(cu) = G [°° — p(v)dv. 

v ; Jo T(v) + iu ' 

This equality together with Eq. (12) implies that the steady-state storage G'(oo) and loss G"(uo) 
shear moduli read 



G 



L 



POO 



o rgexp(-2w) +uj 2 
00 r exp(— v)uj 



exp(— 2v ) + uj'' 



■p(v)dv, 
■p(v)dv. 



(27) 



To fit the experimental data, we adopt the random energy model (Derrida, 1980) with the 
quasi-Gaussian distribution function p{y ), 



p(y) = Po exp 



(v - Vf 



2S 2 



{v > 0), p(v) = (v< 0), 



(28) 



where V and E are adjustable parameters (an apparent average activation energy and an 
apparent standard deviation of activation energies, respectively), and the constant p is found 
from the normalization condition _ 

poo 

/ p( v )dv = 1. (29) 

JO 

The average activation energy for separation of strands V* and the standard deviation of acti- 
vation energies £* are determined by the conventional formulas 



V* 



roo 

I vp(v)dv, 
Jo 



(v - V^) 2 p{v)dv 



(30) 



Governing equations (27) and (28) involve three material constants: (i) the instantaneous shear 
modulus G, (ii) the average activation energy for rearrangement of strands in a network V*, 
and (iii) the standard deviation of activation energies £*. For shear oscillatory tests performed 
at various temperatures T, these quantities become functions of T. Our purpose now is to find 
these parameters by matching the observations for G'{uj) and G"{uj) depicted in Figures 1 to 4. 
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4 Fitting of observations 



Each set of the experimental data is fitted separately. We fix some intervals [0, V max ] and 
[0, S max ], where the "best-fit" parameters V and £ are assumed to be located, and divide these 
intervals into J subintervals by the points = iAV and = jAS = 1, . . . , J — 1) 
with Ay = V max /J and AS = S max /J. For any pair {V {i \ £ w }, the coefficient p in Eq. (28) 
is determined from Eq. (29), where the integral is evaluated numerically by Simpson's method 
with 400 points and the step Av = 0.1. The integrals in Eq. (27) are calculated by using the 
same technique. The shear modulus G is found by the least-squares method from the condition 
of minimum of the function 



where the sum is calculated over all frequencies oo m at which the data were collected, G' and 
G" xp are the storage and loss moduli measured in a test, and G' num and G'^ um are given by Eq. 
(27). The "best-fit" parameters V and £ are determined from the condition of minimum of 
the function F on the set {V^\ S^}. After finding the "best-fit" values V (i) and this 
procedure is repeated twice for the new intervals \V { - i ~ 1 \ and [E^ _1 \ E^ +1 ^], to ensure 

an acceptable accuracy of fitting. Figures 1 to 4 demonstrate good agreement between the 
experimental data and the results of numerical simulation. 

After finding the instantaneous shear modulus G of iPP and LDPE, we calculate the pa- 
rameter l e by using Eq. (5) and plot this quantity versus temperature T in Figure 5. To 
demonstrate the level of discrepancies between the values of l e measured on different speci- 
mens, experimental data are presented for two sets of observations on iPP. The dependence 
/ e (T) is approximated by Eq. (8), where the coefficients (m = 0, 1) are determined by 
the least-squares technique. Figure 5 shows that Eq. (8) provides fair approximation of the 
observations. The average contour length between entanglements l e grows with temperature T. 
The rate of increase in l e {T) is more pronounced for iPP than for LDPE. 

Given adjustable parameters V and E determined at each test temperature T, we calculate 
the dimensionless average activation energy for separation of strands from temporary junctions 
V* and the dimensionless standard deviation of activation energies £* by formulas (30), where 
the integrals are evaluated numerically. The average activation energy V* and the standard 
deviation of activation energies £* are found from the equations similar to Eq. (11), 



The average activation energy V* is depicted versus temperature T in Figure 6. The experi- 
mental data are approximated by the linear equation 



where the coefficients V m (m = 0, 1) are determined by the least-squares method. Figure 6 
demonstrates that Eq. (31) correctly describes the evolution of V* with temperature. The 
average activation energy for rearrangement of strands V* grows with temperature for both 
polymers. The rate of increase is more pronounced for LDPE than for iPP. Given a temperature 
T, the average activation energy for separation of strands in iPP noticeably exceeds that in 




K = k B TV* 




K = (Vo + V 1 T) ■ 10 



19 



(31) 



LDPE. 
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The standard deviation of activation energies is plotted versus temperature T in Figure 
7. The experimental data are approximated by the linear function 



£* = (S + EiT) • 1(T 19 , (32) 

where the coefficients S m (m = 0, 1) are determined by the least-squares technique. Figure 
7 reveals that Eq. (32) adequately describes changes in with temperature. The standard 
deviation of activation energies S* slightly increases with T for iPP and weakly decreases with 
temperature for LDPE. Given a temperature T, the standard deviation of activation energies 
for LDPE substantially exceeds that for iPP. 

To assess the effect of temperature T on the average relaxation time r , we calculate the 
modulus of complex viscosity r\ by the standard formula 



(Gv» 2 +(^r 



(33) 



and find the zero-frequency complex viscosity i] = lim^ vi 00 )- ^ follows from Eqs. (27) and 
(33) that 

G f°° 

T] = — / exp(v)p(v)dv. (34) 
1 o Jo 

The average relaxation time Tq is given by 

Vo 

To= G- 

Substitution of expression (35) into this equality results in 

1 f°° 

r = — / exp(v)p(v)dv. (35) 
1 o Jo 

For any temperature under consideration T, we calculate the integral in Eq. (35) by Simpson's 
method with 400 points and the step Av = 0.1. The distribution function p(v) is given by Eq. 
(28) with the parameters V and S found by fitting the experimental data for G'{uj) and G"{uj). 
The average relaxation time r is plotted versus temperature T in Figure 8. The experimental 
data are approximated by the Arrhenius dependence 

t = r,exp(-^), (36) 

where r* is the average relaxation time at elevated temperatures (T — > oo), E is the activation 
energy and R is the universal gas constant. To match the observations, we present Eq. (36) 
in the form 

(i) 

lnr = r (0) + ^ (37) 

with 

Tf=lnr„ T W = ! (38) 

Figure 8 demonstrates that Eq. (37), where the coefficients Tq 71 ^ (m = 0, 1) are determined by 
the least-squares technique, adequately describes the experimental data for both polymers. 
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5 Discussion 



Figures 5 shows the average contour length between entanglements l e strongly increases with 
temperature for both polyolefins under investigation. The characteristic temperature for disen- 
tanglement T e of LDPE exceeds that of iPP by about twice. To calculate the coefficient a e in 
Eq. (7), we use Eq. (8) and the results of numerical analysis and find a c = 0.023 Kr 1 for iPP 
and a c = 0.013 K _1 for LDPE. This means that l e changes rather weakly with temperature 
for a highly branched polymer melt, and it is strongly affected by temperature for a melt with 
relatively short branches. 

This finding may be explained based on the following scenario. Slightly above its melting 
temperature T m , polypropylene chains are closely packed, which implies that the number of 
macromolecules that intersect a volume occupied by an individual chain is relatively large (a 
high average number of entanglements per chain and a small distance between entanglements 
l e ). With the growth of temperature, thermal oscillations of segments induce loosening of this 
packaging and an increase in the occupied volume. Some macromolecules that crossed the 
occupied volume for an individual chain near the melting point are forced to leave this volume 
with an increase in T. This results in a substantial decrease in the number of entanglements 
per chain and a strong exponential growth of l e . 

On the contrary, due to the presence of long branches, packaging of chains in LDPE at the 
temperatures T close to its melting temperature T m is relatively loose. This poor packaging is 
associated with formation of physical cross-links (knots) between long branches and between 
long branches and the backbones of polyethylene chains (in addition to entanglements between 
backbones). The growth of temperature induces partial disentanglement of junctions formed 
by long branches, which implies that the main chains become more closely packed (the latter is 
reflected by the fact that the coefficient of thermal expansion of LDPE is negative). Although 
the total number of entanglements between chains is reduced with temperature, the rate of 
increase in l c (T) in LDPE is weaker than that in iPP, in agreement with the experimental data 
depicted in Figure 5. 

Figure 6 demonstrates that the average activation energy for separation of strands from 
temporary junctions V* slightly grows with temperature for iPP. This increase may be associated 
with the fact that disentanglement of chains driven by the growth of temperature begins with 
relatively weak junctions (that are characterized by small activation energies). The latter means 
that the average activation energy of an ensemble of chains increases due to disappearance of 
"weak" junctions untangled under heating. 

The rate of growth of the average activation energy V* of LDPE with temperature T is 
substantially stronger than that of iPP. This may be explained by the fact that an increase 
in temperature induces breakage of physical cross-links between long branches and between 
long branches and backbones in LDPE (which are characterized by relatively low activation 
energies). As a result of thermally-induced destruction of "weak" junctions, only relatively 
strong knots between backbones survive. Because these "strong" entanglements have high 
activation energies, the average activation energy of an ensemble pronouncedly increases with 
temperature. 

According to Figure 6, the average activation energy V* of iPP substantially exceeds that 
of LDPE. This conclusion appears to be natural, because the activation energy of a junction 
is determined by local properties of entangled chains. The strength of a knot formed by two 
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chains is relatively low for LDPE (interactions between backbones of linear chains), and it is 
noticeably higher for iPP due to additional interactions between methyl side-groups belonging 
to different chains. 

Figure 7 shows that the standard deviations of activation energies X* of iPP and LDPE 
are weakly affected by temperature. The parameter slightly increases with temperature for 
iPP (due to partial disentanglement of backbones), and it slightly decreases with T for LDPE 
(driven by breakage of weak junctions between long branches). Given a temperature T, the 
standard deviation of activation energies of LDPE substantially exceeds S* of iPP. This finding 
seems quite natural, because the presence of "weak" physical cross-links between long branches 
implies a broad distribution of activation energies for rearrangement of strands. 

Figure 8 reveals that the average relaxation time r grows with temperature for both poly- 
mers under consideration. Given Tq~\ we calculate the activation energy E from Eq. (38) and 
find E = 10.4 kcal/mol for iPP and E = 18.1 kcal/mol for LDPE. These values are in good 
accord with the activation energies provided by other researchers for polypropylene: E = 9.3 
kcal/mol (Eckstein et al, 1998), E = 9.7 - 10.0 kcal/mol (Pearson et al., 1988), E = 10.0 
kcal/mol (Fujiyama et al., 2002), and low-density polyethylene: E — 9.1 — 13.2 kcal/mol 
(Wood-Adams and Costeaux, 2001), E = 13.5 kcal/mol (Qiu and Ediger, 2000). 

6 Concluding remarks 

Two series of torsional oscillation tests have been performed on melts of isotactic polypropylene 
and low-density polyethylene in the range of temperatures between 170 and 250 °C (iPP) and 
between 120 and 190 °C (LDPE). 

With reference to the concept of transient networks, constitutive equations have been de- 
veloped for the viscoelastic response of polymer melts at isothermal three-dimensional de- 
formations with small strains. The melt is treated as an equivalent transient network of 
strands bridged by temporary junctions. Its time-dependent behavior is modelled as thermally- 
activated separation of active strands from their junctions and attachment of dangling strands 
to the network. Stress-strain relations for an equivalent heterogeneous network of strands 
(where different junctions have different activation energies for rearrangement of strands) have 
been derived by using the laws of thermodynamics. These equations involve three material 
parameters that are determined by matching the experimental data for the storage and loss 
moduli as functions of frequency of oscillations. Fair agreement is demonstrated between the 
observations and the results of numerical simulation. 

The following conclusions are drawn: 

1. The average contour length between entanglements l c grows with temperature T for both 
polyolefins under consideration. Thermally-induced changes in l c (T) are correctly de- 
scribed by Eq. (7). The relative rate of increase in l e (T) is higher for iPP than for 
LDPE. 

2. The average activation energy for rearrangement of strands V* grows with temperature 
due to breakage of "weak" junctions between chains. The rate of increase in V* for 
LDPE exceeds than for iPP (due to thermally-induced destruction of knots between long 
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branches). Given a temperature T, the average activation energy of iPP is substantially 
higher than that of LDPE (due to local interactions between methyl side-groups). 



3. The standard deviation of activation energies for separation of active strands E* is weakly 
affected by temperature. The value of for LDPE noticeably exceeds that for iPP (which 
reflects the fact that the presence of "weak" physical cross-links between long branches 
in polyethylene results in a very broad distribution of activation energies). 

An explicit expression (35) is derived for the average relaxation time r . It is demonstrated 
that the activation energies of iPP and LDPE calculated by using Eqs. (35) and (36) are close 
to those reported in the literature. 
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List of figures 



Figure 1: The storage modulus G' versus frequency to. Circles: experimental data on iPP at the 
temperatures T = 190, 210, 230 and 250 °C, from top to bottom, respectively. Solid lines: results on 
numerical simulation 

Figure 2: The loss modulus G" versus frequency u. Circles: experimental data on iPP at the 
temperatures T = 190, 210, 230 and 250 °C, from top to bottom, respectively. Solid lines: results on 
numerical simulation 

Figure 3: The storage modulus G 1 versus frequency to. Circles: experimental data on LDPE at the 
temperatures T = 120, 130, 140, 150, 160, 170, 180 and 190 °C, from top to bottom, respectively. 
Solid lines: results on numerical simulation 

Figure 4: The loss modulus G" versus frequency lo. Circles: experimental data on LDPE at the 
temperatures T = 120, 130, 140, 150, 160, 170, 180 and 190 °C, from top to bottom, respectively. 
Solid lines: results on numerical simulation 

Figure 5: The parameter l e versus temperature T. Symbols: treatment of observations on iPP 
(unfilled circles) and LDPE (filled circles). Solid lines: approximation of the experimental data by 
Eq. (8). Curve 1: = 1.09, = 1.04 • 10~ 5 , T = 43.2. Curve 2: = -0.09, = 1.30 • 10~ 3 , 
T = 77.2 

Figure 6: The average activation energy T4 versus temperature T. Symbols: treatment of observa- 
tions on iPP (unfilled circles) and LDPE (filled circles). Solid lines: approximation of the experimental 
data by Eq. (31). Curve 1: V = 4.42 • 10 3 , V x = 16.41. Curve 2: V = -1.06 • 10 4 , V\ = 40.22 

Figure 7: The standard deviation of activation energies S* versus temperature T. Symbols: treat- 
ment of observations on iPP (unfilled circles) and LDPE (filled circles). Solid lines: approximation of 
the experimental data by Eq. (32). Curve 1: S = 1.01 • 10 3 , Si = 1.79. Curve 2: S = 3.36 • 10 3 , 
Si = -9.52 

Figure 8: The average relaxation time To versus temperature T. Symbols: treatment of observations 
on iPP (unfilled circles) and LDPE (filled circles). Solid lines: approximation of the experimental data 
by Eq. (37). Curve 1: r = -13.74, n = 5.24 • 10 3 . Curve 2: r = -22.18, n = 9.12 • 10 3 
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